Can substructure in the Galactic Halo explain the ATIC and PAMELA results? 
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Recently, ATIC and PAMELA measured an anomalously large flux of leptonic cosmic rays which may arise 
from dark matter self-annihilation. While the annihilation signal predicted for a smooth halo is 10^ — 10"' times 
smaller than the measured excess, the signal can be boosted by the presence of subhalos. We investigate the 
feasibility of large boost factors using a new Monte Carlo calculation technique that is constrained by previous 
simulation work on halo substructure. The model accounts for the observed decrease in the amount of substruc- 
ture with decreasing halo mass and the scatter in halo structural parameters such as the density concentration 
parameter. Our results suggest that boost factors of ~ 10^ are ruled out at > 14(t. We conclude that substructure 
alone, at least with commonly assumed annihilation cross-sections, cannot explain the anomalous flux measured 
by ATIC and PAMELA. 
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Recent measurements by ATIC H and PAMELA lH offer 
the tantalizing prospect that dark matter has been discovered, 
albeit indirectly. Both experiments have reported an anoma- 
lously large flux in leptonic cosmic rays at energies above 
10 GeV. Dark matter candidates, such as WIMPs (Weakly 
Interacting Massive Particles which arise in theories of su- 
persymmetry) or Kaluza-Klein particles, can annihilate and 
produce cosmic rays {sS]. The excess flux might comprise sec- 
ondary particles from annihilation events in the Galactic halo, 
a possibility that is now attracting considerable attention. 

Though the observed anomalous flux is 100-1000 times 
larger than the flux predicted from a smooth Milky Way halo 
with a standard thermal WIMP particle, the presence of sub- 
halos can boost the signal since the annihilation rate is pro- 
portional to the square of the density. But while there is little 
doubt that dark matter halos are clumpy, the actual boost fac- 
tor is a matter of some debate. In this paper, we take a critical 
look at the underlying assumptions in boost factor calculations 
and discuss the implications of our analysis for the ATIC and 
PAMELA results as well as current and future observations by 
the Fermi Gamma-Ray Large Area Space Telescope (GLAST) 

Halos in cosmological N-body simulations host numerous 
subhalos which in turn host their own subhalos |@]. The 
distribution of subhalos can be summarized by the mass func- 
tion, dN/dlnf, where / = m/Mh is the relative mass frac- 
tion of a subhalo of mass to in a host of mass Alh - Note that in 
this paper, a host can refer to either a halo or a subhalo. The 
mass function appears to be well-characterized by a power- 
law ill[i3, 

dN/d\nf = Af~^, (1) 

for / < 10^^. While current simulations probe dA^/d In / 
for / > 10"^, a Galactic-mass halo with neutraUno dark mat- 
ter should have subhalos down to f ^ 10~^* ifTTIl . with the 
possibility of up to nine nested levels of substructure. 

The luminosity of annihilation secondaries (electrons, 
positrons, or photons) from a halo or subhalo of mass 



can be written in the form 

where Ve± are the branching ratios for electron-positron or 
photon secondaries, a is the total annihilation cross section, 
and V and are the relative speed and mass of the annihi- 
lating particles. The quantity C{Mh) is the volume integral of 
where p is the dark matter density. The boost factor due to 
subhalos is defined implicitly by the relation 

C{Mh) = [l + B{Mh)]C{Mh). (3) 

Here C{Mh ) represents the volume integral for a smooth halo, 
given by the volume integral of the smooth density field fP. 
The boost factor is calculated recursively through the expres- 
sion 

BiM,)^jl^J-^il + Bin.)nrn)dlnf, (4) 

where m = fMh- 

Gamma-ray secondaries travel directly from the source to 
the observer and therefore the flux measured at Earth due 
to a source located at a distance d is L^/Aird^. In contrast, 
charged particles are scattered by the Galactic magnetic field 
and therefore the arrival direction of the secondaries is essen- 
tially independent of the direction to the source. Furthermore, 
charged particles lose energy en route from the source to the 
observer It is therefore useful to define an effective (energy- 
dependent) total boost due to all subhalos (see Ref. lfl2ll for 
details): 

B,n{E;Mh,E,,^Q) = {I - ftf 
X0)esub(A4) 

Here ft is the total mass fraction in subhalos and again m = 
fMh- The first term accounts for the fact that p is defined 
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FIG. 1: Compilation of a (top) and /42 (bottom) from the studies indicated in right-hand key, where the year of the study decreases as one 
goes down the column. For both the upper and lower left panels, marker colors indicate different studies while the marker type indicate the 
algorithm used to identify subhalos. Studies which examined subhalo hosts instead of field halo hosts are indicated by internal black filled 
star. Colored horizontal error bars indicate mass range examined in study and black dashed vertical lines going from left to right indicate mass 
dwarf spheroidal galaxies, galaxies, and galaxy clusters. In the upper and lower right panels, we plot our phenomenological models of a and 
/42 with the thick and thin lines for model I and respectively. For a, the mean, upper and lower values are given by the solid, dashed and 
dotted lines respectively. For (42, the mean of the lognormal distribution is given by the solid line and +2a by the dashed, —2a by the dotted 
lines. Top panels: vertical black and grey eiTor bars indicate uncertainty in the published and fit-by-eye values respectively. Bottom panels: 
black vertical error bars indicate variation between similar mass hosts. 



such that J pdV = Mh, that is the halo is smooth, but only 
(1 — /t) of the host's mass is in the smooth radial component. 
Essentially (1 — /t)^ is a normalization term. 

In Eq. (|5]l, the effect of subhalos has been broken up into 
two terms by assuming that the volume distribution of subha- 
los is independent of mass distribution. The first term in the 
numerator, ^sub, accounts for the propagation of cosmic rays 
originating from subhalos and is given by 

e,ub(S;£;^,X0) y"G(£;,X0;£;„x,)^d'x,, (6) 

where the host term has been normalized by the solar den- 
sity pq. Here the propagation of the secondaries from 
an initial position x, with initial energy Ei to the Earth 
with final energy E is explicitly accounted for using the 
Green function G{E,x.q, Ei,Xi). The propagation is a 
diffusive process and consequently G{E,Xq, Ei,Xi) oc 
exp [— (xi — X0)^/A|)(i?; £'i)] where is the diffusion 



length. For electron-positron pairs produced by a 100 — 
1000 GeV WIMP, Xd few kpc for E > 100 GeV and de- 
creases monotonically as E Ei (see Ref. |fl2tl for further 
discussion). In the definition of Qsuh, we are effectively treat- 
ing subhalos as point sources, which is a reasonable assump- 
tion as w 90% of the flux originates from the central region 
of a (sub)halo. As we are interested in subhalos that are nu- 
merous enough to enhance the diffusive background, we are 
generally concerned with subhalos with masses of < 10* Mq. 
These halos have central regions that are < kpc in radius. We 
also assume in Eq. ^ that subhalos trace the host's dark mat- 
ter distribution. 

The second term in the numerator, ^sub, is the total contri- 
bution of subhalos to the annihilation flux and is given by 

UhiMh)^ J ■^C{m)[l + B{m)]dlnf. (7) 

This is equivalent to the cosmic ray source term. The main 



3 



difference between our work and previous studies 
is that we examine the enhancement that arises explicitly due 
to the entire subhalo hierarchy, that is subhalos, subsubhalos, 
etc., by incorporating the boost factor. 

The ^haio term in denominator of Eq. (|5]l is the host's con- 
tribution to the cosmic ray flux which also accounts for the 
propagation of cosmic rays and is given by 

£.h^io{E;E„XQ) ^ J G{E,XQ;E„x,)p^{x,)d^:>c,, (8) 

Here again the propagation of the secondaries is explicitly ac- 
counted for using the Green function G{E, Xq, Ei, x^). 

The effective boost in cosmic ray flux can be greatly en- 
hanced by the boost factor. The boost factor depends sensi- 
tively on the subhalo mass function. Using the often quoted 
values a 1 and A 0.033 iB, one finds B ^ 30 for 
7-rays with a Galactic-mass halo. The key assumption in ob- 
taining this result is that a and A are independent of the host 
mass and apply to all scales and levels in the subhalo hierar- 
chy. Our goal is to examine the validity of this assumption, 
in short, to test whether a and A depend on the mass of the 
host. To do so, we compare estimates for a and A over a wide 
range of published simulations. 

Though most of the results in our study are based on sim- 
ulations which assume a standard ACDM cosmology, results 
from our own simulations of structure formation in scale-free 
cosmologies lfl6ll as well as results from simulations of warm 
dark matter cosmologies UtIi are also included. To assign 
an effective mass to hosts in these non-standard simulations 
we appeal to the halo formation process. In the hierarchi- 
cal clustering paradigm, halos arise from primordial density 
fluctuations where small-scale structures form first and then 
merge together to form larger and larger objects. This process 
is governed by the power spectrum of the perturbations, P{k), 
where k is the comoving wavenumber of a perturbation mode. 
Alternatively, one can describe the perturbation spectrum by 
the mass variance, 

-^(M) = J ^Pik)WHk,M), (9) 

where W{k, M) is the window function enclosing a mass 
M oc k~^ within a comoving volume of radius r = 27r/fc. 
This quantity is simply the average rms overdensity of a 
sphere enclosing a mass M. 

In scale-free cosmological models the primordial power 
spectrum is a power-law function of wavenumber k, that is, 
P oc fc" where n is referred to as the spectral index. This 
form of the power spectrum leads to a mass variance obeying 
(ihiCT^/(ilni\jf = — (n + 3)/3. Though the power-spectrum 
in a ACDM cosmology is more complicated, we can define an 
effective spectral index, nsfi{M) = —3{d\na^ / d\n M + 1). 
Thus, the spectral index in a scale-free cosmology may be 
used as a proxy for the (ACDM) mass scale [160 by set- 
ting n = neff{M) and solving for M. In this work we take 
= 0.25, Qa = 0.75, h = 0.73, = 0.9 and = 1 as 
our reference cosmology. 

The use of different algorithms to identify subhalos intro- 
duces a further complication. These algorithms range from 



SUBFIND ifisll . which associates subhalos with local density 
peaks, to 6DFOF d [H, which searches for clustering in 
6D phase space. In addition, some (but not all) researchers 
apply an unbinding criterion which removes unbound parti- 
cles when estimating the mass of a subhalo. One reason for 
such a variety of methods is that subhalo identification is an 
ill-defined problem. 

To characterize the amount of substructure, we introduce 
/42 in place of A, where 

-2 In 10 

U.^ J f^di^f, (10) 

-41n 10 

is the fraction of mass in subhalos with 10~^ < / < 10^^. 
In Fig. [T] we plot estimates of a and /42 from a large sam- 
ple of studies. When available, we plot published values of 
a. Otherwise, we fit-by-eye for the mean a and estimate the 
uncertainty. The mean a is then used to determine /42. Error 
bars for /42 are shown when a study reports the variation in 
the subhalo number or mass fraction across several hosts of 
the same mass. 

The figure reveals a large amount of scatter in a and /42. 
For galactic to cluster masses, estimates of a range from 0.7 
and 1.1, which is generally within the estimated uncertainty 
in a for a single measurement and the variation in a with red- 
shift for an individual halo lfl6l I20I1 . The scatter in log /42 
is roughly constant with Mh with most points lying within 
0.15 dex of the mean. The scatter may be due to variations 
in the mass accretion histories of the halos. The halo's mass 
accretion history is important not only for defining the density 
parameter ll2lll but it is also responsible for delivering new 
substructure within a halo. 

Figure[T]also suggests that a and /42 decrease with decreas- 
ing Mh. The observed A//i -dependence should be treated with 
some caution. For example, a linear fit to a{Mh) for Mh be- 
tween galaxy and cluster masses is consistent at the Icr-level 
with a = constant. Furthermore, at smaller host halo masses 
(neff closer to —3) the subhalo mass function is more sensi- 
tive to the type of subhalo-finding algorithm used lfl6ll and in 
particular, whether a binding criterion is applied. For smaller 
Mh, subhalos tend to be less gravitationally bound. Studies 
that do not correct for unbound particles, such as lfl9ll . may 
overestimate subhalo masses and therefore overestimate a and 
/42. Even for larger Mh, the application of a correction for 
unbound particles decreases a by 0.1 — 0.2 (see, for exam- 
ple, Ref. 1I24I23II). It should also be noted that the data points 
from Ref. ll24ll at subgalactic masses are for subhalo hosts. 
This study found that the mass fraction of subsubhalos in sub- 
halos tends to be smaller than the mass fraction of subhalos 
in halos. It may well be that a depends on the host's level in 
the subhalo hierarchy, that is a of a field halo may differ from 
that of a similar mass host that sits within the virial radius of 
a larger halo. 

Most boost factor calculations assume a and /42 are inde- 
pendent of mass, that is, a{Mh) ~ ao, f42{Mh) = fo- How- 
ever, this assumption does not capture the trends in Fig.[Tl es- 
pecially at lower Mh, and we therefore propose the following 
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TABLE I: Summary of (a, fi2) model parameters 



Model # 


Fig.|2]color 


OLo 




/o 


u 


0-u 


g 


1.0 





0.04 








r 


0.9 





0.04 





0-1 


b 


0.7 





0.04 





1-u 


g 


1.0 


0.046 


0.034 


10 


1 


r 


0.9 


0.12 


0.034 


10 


1-1 


b 


0.7 


0.46 


0.034 


10 



phenomenological model: 

a{Mh) = ao + a„ log[neff(A4) + 3], (11) 
\ogU2{Mh) = log/o + log/„ log[neff(M,,) + 3]. (12) 

We account for the scatter in a by considering three values for 
each Qfo and a„ - the mean and the upper and lower envelop 
values as given in Table U For f^^ we fit the mean by eye and 
assume that it follows a lognormal distribution. This choice 
is motivated by the fact that other bulk halo parameters, such 
as the concentration parameter, follow lognormal distributions 

We calculate the boost factor using Eq. (01, under the as- 
sumption that the subhalo mass function is a power-law with 
an index a ( Mh ) ■ The amplitude is normalized using /42 ( Mh ) 
and a{Mh) subject to the condition that the total fraction of 
mass in subhalos is less than one. The mass function ex- 
tends down to a mass mo, the minimum CDM halo mass, 
which depends on the properties of the dark matter particle. 
Typical values for TOq when dark matter is a neutralino in 
the Constrained Minimal Supersymmetric Standard Model are 
10"^ - 10~*^ Mq UH], though the minimum CDM halo mass 
can vary between 10~^^ — 10^'' Mfrj26t]. We use the more 
probable value of ttIq = 10^^ Mq jllll . We also assume that 
a host cannot contain subhalos with / > 10^^ as suggested 
by most studies. As a consequence, the smallest host mass 
that contains subhalos is lOOmo. Finally, we only go to three 
levels in the subhalo hierarchy since deeper levels change the 
boost factor by < 1%. As a general observation the mass flux 
tends to come from the first two levels of the subhalo hierar- 
chy. 

Due to the condition that ft < 1, the boost factors begin 
to saturate for large /42 and small mo- The limiting values of 
/42 for a given a and mo/Mh is given in Table HIl This table 
shows that the limits are only important for a > 1 and only 
when mo/Mh is large, for example, when jtIq = 10^^^ Mq 
for Mh = 10^2 Mq. For our choice of jtIq = lO^*^ Mq, 
galactic mass hosts, and our subhalo mass function models, 
this saturation occurs at a /42 that is still « 3cr away from the 
mean even at a = 1.0. 

The volume integral of in Eq. Q is given by 

C{Mh) = ^MhC^Mh)gic), (13) 

where py is the virial density of a halo, c = ry/rg is the con- 
centration parameter, is the profile's effective radius, and 
the dimensionless function g(c) depends on the form of the 



TABLE IL Limits on /42 due to /t < 1 
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density profile. We consider the commonly used iNFWl ll27ll 
profile, which is described by two parameters, the characteris- 
tic radius and density ps. We also examine the Einasto pro- 
file, which appears to be a better fit to halos from cosmological 
simulations fl% . The Einasto profile has an extra parameter 
describing the radial logarithmic slope of the density profile 
which appears to have a mass dependence ll28ll . We find that 
an Einasto profile increases the boost factor for subgalactic 
masses by < 50%. To determine the concentration parame- 
ter we use the two-parameter model presented in Ref. Izsll for 
c{M) as this model provides a better match to simulation data 
than the often-used model of Ref. ll29l] . Using the Ref. ll25ll 
model for c{M) instead of the Ref. 112911 model reduces B by 
< 10%. 

We calculate 10^ random realizations of B for a given Mh 
where c and /42 are sampled from lognormal distributions 
with criog/42 = 0.15 dex and ciogc — 0.10 dex ll30ll . Due to 
the recursive nature of Eq. (HI, variation in c not only affects 
C of the host but that of the host's subhalos. The distribution 
of B{Mh) across the 10"* realizations appears to be well char- 
acterized by a lognormal distribution. We fit this distribution 
for the mean and the variance. 

The results are plotted in Fig.|2] We particularly highlight 
the mass scales of the Galaxy's dwarf spheroidal satellites as 
these objects may be the best candidates for searches of 7-rays 
secondaries |5]. The left panel shows how a few key parame- 
ters individually affect B{Mh)- As a reference model, we use 
mo = 10"^ Mq, a = 1.0, /42 = 0.067, and anjNFW| density 
profile with the Bullock prescription for c(M) II15I1 . Note that 
in Ref. ifTsll . the boost factor is calculated for subhalos below 
the mass resolution of their Via Lactea II simulation, which 
is ^ 10^ Mq. This amounts to calculating Eq. (HI with the 
upper limit given by 10^ Mq/M/j. Decreasing a from 1.0 to 
0.9, or using the mean values from our new phenomenological 
model for /42 decreases -B(A/dSph) by ~ 4. The introduction 
of scatter in c and /42 increases the mean boost factor slightly, 
though the amount depends on the exact form and width of the 
distributions. For our choices, -B(-A/dSph) increases by w 30%. 
Decreasing mo to 10^^ Mq also increases the boost factor by 
K, 30% while using an Einasto profile increases the boost fac- 
tor by < 20% in this case. 

The other two panels of Fig. |2] show the peak of the distri- 
bution and the 2a contours from the models listed in Table U 
We find increasing a increases the boost factor by w 2 for 
dSph galaxies, though the differences are generally within the 
variation caused by c and /42. For model 0, the mean a gives 
-B(MdSph) ^ 0.6lJ;4, though for a = 1.0, B - 8 is within 
the 2(7 envelop. Model 1 reduces the boost factors such that 
i3(AidSph) ^ 2 lie outside the 2cr envelop with the mean a 
giving B(AfdSph) ~ 0.2. 
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10^ 10' 10° 10' 10'° 10" 10'^ 10° 10' 10° 10' 10'° 10" 10'^ 10° 10' 10° 10° 10" 10" 10'^ 

FIG. 2: Boost factors. Left panel shows how B{Mh) changes when a single parameter is changed. Middle (right) panel shows model 0(1) 
with thick red, green and blue lines indicating mean boost factor using the mean, upper, and lower values of a respectively. Colored hatched 
contours indicate the 2a region arising from the variation in lege & log/42. The gray hatched region outlines the mass range for dwarf 
spheroidal satellite galaxies of MW. 



Thus, subhalos in the Galaxy's satellites are unlikely to 
greatly enhance the 7-ray flux. Previous estimates of the num- 
ber of satellites that could be detected via their 7-ray flux by 
GLAST, such as those by Ref. Q [ist], are probably overly 
optimistic. Even our revised calculations may overestimate 
the boost factor as it appears that subhalos have smaller 
than similar mass halos. Consequently, satellites might even 
have their boost factors reduced by a factor of > 2. So far, 
GLAST has detected numerous 7-ray sources, none of which 
are convincing dark matter annihilation signals llsill . 

We now examine the consequences of such boost factors for 
cosmic rays. To determine the energy dependence in Eq. dUl , 
we use the mean propagation parameters listed in Ref. 111211 
and assume that the dark matter particle annihilates directly 
to monoenergetic e^. For a 700 GeV thermal dark matter 
particle with {av) = 3 x 10^^^ cm^s~^, the diffusion length 
is 1.3 kpc at i5 = 300 GeV and monotonically decreases to 
0.1 kpc at E = 69 GeV. We also assume that the Milky Way 
halo has an INFWI density profile with a characteristic radius 
Ts = 20 kpc and a solar density of pq = 0.43 GeV cm~'^. 

Recall that in Eq. (|5]l we assumed that the subhalo volume 
distribution is the same as the halo's density profile. How- 
ever, simulations appear to show that subhalos do not trace the 
host halo's radial density profile ll24ll35ll . Instead tidal strip- 
ping destroys subhalos in the center, resulting in a radially 
anti-biased distribution, at least for subhalos with masses of 
> 10^ Mq. Very few subhalos are thus found within ~ 20 kpc 
of the halo center. This implies is that subhalos are unlikely 
to contribute to the observed flux as they are too far away. We 
argue that it is unreasonable to assume that this anti-biased 
radial distribution continues down to very small masses since 
these subhalos should become increasingly less susceptible to 
tidal disruption as their mass decreases and therefore should 
be able to survive to smaller radii. It is also worth noting 
that it is difficult to identify substructure in the central, high 



density regions of halos with current group finder algorithms. 
Furthermore, Ref. Issll found that though tidal fields affect 
the mass of subhalos, their annihilation luminosity is less af- 
fected and more closely follows the host halo's radial density 
profile. Thus, we compromise and only include contributions 
from subhalos with m < 10^ M©. Considering only low 
mass subhalo hosts reduces the contribution of substructure 
by w 60% for both models and 1 . This may still be an op- 
timistic assumption since interactions with baryonic matter, 
that is the Galactic disk and stars, might further reduce the 
inner subhalo abundance. 

In general we find B^ff{E > 300 GeV) ^ 1. Even the most 
optimistic model, 0-u, gives Beff(-E > 300 GeV) = l-ilol- 
Here the variation in Bgff is due primarily to the variation 
in /42 for the Galactic halo based on our phenomenological 
model, though we also account for the variation in the boost 
factors and concentration parameters of subhalos. This model 
excludes the i?eff ~ 200 required to explain the ATIC signal 
at the ~ 14(T level. These results are in agreement with sim- 
ilar studies lfl2[ [13I [14I1 . We have improved on earlier results 
by explicitly including several levels in the subhalo hierarchy 
that earlier work neglected. 

The High Energy Stereoscopic System (HESS) ll34ll and 
Fermi GLAST |01 found different though not necessarily con- 
tradictory results. Both instruments detect an excess, though 
not as large as that observed by ATIC, and neither instrument 
observes the pronounced peak at ~ 700 GeV seen by ATIC. 
The Fermi GLAST measurements require a boost of ^ 150 at 
E = 300 which monotonically decreases to « 70 at = 600. 
This is still rejected by our model at the 10 — 12ct level. 

Some studies have suggested that it is possible to repro- 
duce the observed flux, particularly the ATIC biimp in the en- 
ergy spectrum, with a single nearby subhalo ll4 I32I1 . The 
required annihilation rate is ~ 10"^^ s^^ for a thermal WIMP 
with 771^ ~ 100 — 1000 GeV. Neglecting the boost factor from 
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deeper levels in the subhalo hierarchy, a subhalo with a mass 
of > 10^ Mq at ~ 1 kpc is required. However, Ref. ll33ll find 
based on a numerical simulation that the probability of such 
a large subhalo within ~ 8 kpc of the Galactic center is ex- 
ceedingly low, p < 10^^. This clearly appears to be highly 
unlikely. Incorporating the boost factors, and thereby incor- 
porating the entire subhalo hierarchy, does not drastically re- 
duce the minimum mass required to explain the ATIC peak 
and thereby the probability of a nearby clump. 

The anomalous cosmic ray flux observed by ATIC and 
PAMELA have sparked a flurry of interest as there is the pos- 
sibility that these instruments may have indirectly detected an- 
nihilating dark matter However, the observed amplitude can- 
not be explain with standard thermal dark matter. A number of 
studies have invoked subhalos to explain the large amplitude 
of this flux. Our work effectively ends this line of thought. 



We find that subhalos alone are unlikely to account for the 
anomaly and are strongly ruled out at the 14a level. This still 
leaves the possibility that the flux is due to instrument errors, 
astrophysical phenomena, or, more intriguingly, perhaps ex- 
otic theories of dark matter are required 1361, ISTT . 
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